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We have developed a relativistic coupled-cluster theory to incorporate nuclear spin-dependent 
interaction Hamiltonians perturbatively. This theory is ideal to calculate parity violating nuclear 
spin-dependent electric dipole transition amplitudes, BlpNCi of heavy atoms. Experimental ob- 
servation of which is a clear signature of nuclear anapole moment, the dominant source of nuclear 
spin-dependent parity violation in atoms and ions. We apply the theory to calculate _E1pnc of Cs, 
which to date has provided the best atomic parity violation measurements. We also calculate -Elp^c 
of Ba + and Ra + , candidates of ongoing and proposed experiments. 
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The effects of parity nonconservation (PNC) in atoms 
occur in two forms, nuclear spin-independent (NSI) and 
nuclear spin-dependent (NSD). The former is well stud- 
ied and experimentally observed in several atoms. The 
signature of the later (NSD) has been observed only in 
one experiment with Cs [Tj and the same experiment has 
provided the most accurate results on NSI atomic PNC 
as well. In an atom or ion the most dominant source 
of NSD-PNC is the nuclear anapole moment (NAM), a 
parity odd nuclear electromagnetic moment. It was first 
suggested by Zeldovich [2] and arises from parity violat- 
ing phenomena within the nucleus. 

One major hurdle to a clear and unambiguous ob- 
servation of NAM is the large NSI signal, which over- 
whelms the NSD signal. However, proposed experi- 
ments with single Ba + ion [3] could probe PNC in the 
Sj/2 — d 5 /2 transition, where the NSI component is zero. 
This could then provide an unambiguous observation of 
NSD-PNC and NAM in particular. The ongoing exper- 
iments with atomic Ytterbium [3] is another possibility, 



To date, the use of PRCC in atomic PNC is limited 
to NSI-PNC. In this letter we report the PRCC theory 
to calculate NSD-PNC in atoms. Such a development is 
timely as the recent experimental proposals on Ba + and 
Ra + |19| and observation of large enhancement in atomic 
Yb [4] shall require precision atomic theory to examine 
the systematics and interpret the results. It must perhaps 
be mentioned that, in an earlier work we had developed 
and calculated electric dipole moment of atomic Hg [20] 
using PRCC theory. 

RCC theory. — In the RCC method, the atomic state 
is expressed in terms of T and S, the closed-shell and 
one-valence cluster operators respectively, as 
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where is the one- valence Dirac-Fock reference state. 
It is obtained by adding an electron to the closed-shell 
reference state, |<J>„) = a||$o)- In the coupled-cluster 
singles doubles (CCSD) approximation = t[ 0) +T 2 (0) 



the 6s 2 % - 6s5d A D 2 transition, to observe NSD-PNC an d = S{ 0) + S^' . The open-shell cluster operators 
with minimal mixture from the NSI component. One cru- are solutions of the nonlinear equations [21] 
cial input, which is also the source of large uncertainty, 
to extract the value of NAM is the input from atomic 
theory calculations. Considering this, it is important to 
employ reliable and accurate many-body theory in the 
atomic theory calculations. 

The coupled-cluster (CC) theory^ E] is one of the 
most reliable many-body theory to incorporate electron 
correlation in atomic calculations. It has been used with 
great success in nuclear [7] , atomic [8HT0] , molecular [TT] 
and condensed matter [T2] physics. In atomic physics, the 
relativistic coupled-cluster (RCC) theory has been used 
extensively in atomic properties calculations, for exam- 
ple, hyperfine structure constants [TO] H3] and electro- 
magnetic transition properties [HI [TS] . In atomic PNC 
calculations too, RCC is the preferred theory and several 
groups have used it to calculate NSI-PNC of atoms [TB1 - 
[T5] . However, the calculations in Ref. [TB] are entirely 
based on RCC with a variation we refer to as perturbed 
RCC (PRCC), where as the calculations in Ref. ^7\M\ 
are based on sum over states with CC wave functions. 



(^\H N + {H^S^}\^ v )=E^ t (^\S[ a) \9 v ), (2a) 
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where £Zn = e _T ' } H-^e Ti ' is the similarity transformed 
Hamiltonian and the normal ordered atomic Hamilto- 
nian H N = H - ($ |#|$o>- And, Ef* = E v - E , is 
the attachment energy of the valence electron. The 
are solutions of a similar set of equations, however, with 
gl ) = 0. A similar set of equations may be derived in the 
case of two-valence systems and use it in the wave func- 
tion and properties calculations of atoms like Yb [22] . 

Perturbed RCC theory. — The perturbed RCC method 
[231 124j , unlike the standard time- independent perturba- 
tion theory, implicitly accounts for all the possible inter- 
mediate states in properties calculations. Consider the 
NSD-PNC interaction Hamiltonian 
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as the perturbation. Here, fi' w is the weak nuclear mo- 
ment of the nucleus and pn if) is the nuclear density. The 
total atomic Hamiltonian is 
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where A is the perturbation parameter. Mixed parity hy- 
perfine states are then the eigen states of Ha- To 
calculate from RCC, we define a new set of clus- 
ter operators T^, which unlike connects the ref- 
erence state to opposite parity states. This is the re- 
sult of incorporating one order of ffp^c an d f° r this rea- 
son we refer to as the perturbed cluster operators. 
Although hyperfme states are natural to i?p|jj=!, clus- 
ter operator TA 1 ) is defined to operate only in the elec- 
tronic space and is a rank one operator. For this define 
H Scc D = ( g f^' H /)/(v / 2) J2i "*PN(r), which operates only 
in the electronic space, so that Hp^c = H ^ C D ■ I. The 
closed-shell exponential operator in PRCC is e T< )+at( 1 
and the atomic state is 



1 + AT (1) • I 
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Similarly, the mixed parity state from one- valence PRCC 
theory is 
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As is one particle and rank one operator, in terms 
of c-tensors 
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where Cj are c-tensor operators. Similarly, the tensor 
structure of Tg is 
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where indicates the two c-tensor operators cou- 

ple to a rank one tensor operator. Based on the tensor 
structures, the perturbed cluster operators are diagram- 
matically represented as shown in Fig. [l] For the doubles 
T2 , to indicate the multipole structure, an additional 
line is added to the interaction line. The cluster opera- 
tors are solutions of the equations 
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FIG. 1. Diagrammatic representation of single and double 
excitation perturbed cluster operators. The short line on the 



interaction line of and S^' is to indicate the multipole 



s(i) 



structure of these operators. 



Where we have used the relations ($p|T (1 ) |$„) = 0, 
and (<FP|TW,S'|<f>„) = 0, as the bra state is valence ex- 
cited. An approximate form of Eq. ([9]), but which con- 
tains all the important many-body effects, are the lin- 
earized cluster equations. This is obtained by consider- 
ing JSnTW 



BnTW, and H^d 
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clcc ^ clcc ' clcc 

We refer to this as the linear approximation and use it 
extensively to check the results. 

^pnc calculations. — If \^f v ) and \^ w ) are atomic states 
of same parity, then the Hp^~, induced electric dipole 
transition amplitude i?lp^° = (* w ||D|^ , l) ), where D is 
the dipole operator. Similarly, the transition amplitude 
within the electronic sector is 
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where D = e T ^De J , is the dressed electric dipole op- 
erator. It is evident that D is a non-terminating series 
of the closed-shell cluster operators. It is non-trivial to 
incorporate T to all orders in numerical computations. 
For this reason D approximated as D w D + DT^ + 
r(0)t D + y(o)t DT (o)_ xhis ca pt urcs a U the important 
contributions arising from the core-polarization and pair- 
correlation effects. Terms not included in this approxi- 
mation are third and higher order in . The expression 
used in our calculations is then 

£bi s c D « ($„|DT« +t(°) , DTW +T( 1 > f DT(°' 
+T d) t D + DT (i) S (o) + T^S^B 
+5 (0)t DT (i) + T (i)t D5 ( 0) + Dg (i) + g(l)t D 

+^ ) t DS( 1 )+S« t D^°)|$ t ,). (II) 

From our previous study of properties calculations |21) , 
we conclude that the contributions from the higher order 
are negligible. 

Coupling with nuclear spin. —To couple Sl^SD with 
nuclear spin I and obtain _Elp^, consider the exchange 
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diagram in Fig. 2^a). It arises from the term T^ ^ DS 2 
in the PRCC expression of £^ C S C D . 

To demonstrate the non-trivial angular integration, in 
hyperfme atomic states, the angular momentum diagram 
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FIG. 2. Examples of .Glp^c diagrams (a) one of the exchange 
diagrams in electronic sector and (b) angular momentum di- 
agram in terms of hyperfine states and the portion within the 
dash lines is the electronic component. 



of the same diagram is shown in Fig. |2jb) . Conventions 
of phase and angular momentum lines of Lindgren and 
Morrison |25j are used while drawing the diagram. The 
portion of the diagram within the rectangle in dashcd-linc 
is the angular momentum part of the electronic sector. 
The evaluation of the angular integral of the electronic 
sector, following Wigner-Eckcrt theorem, is equivalent to 

0„m„|Yjjlf> , DS< 1 >} lj,m„) = (-l)J--"'- 

><i:(-i i *:)«.i{jf , "f}'u.).(u) 



TABLE I. Reduced matrix element, SI^nc, of the 6 2 S 1/2 -> 
7 2 S 1/2 , 6 2 S 1/2 ->■ 5 2 D 3/2 and 7 2 S l/2 6 2 D 3/2 transi- 
tions between different hyperfine states in Cs, Ba + and Ra + 
respectively. The values listed are in units of ieao x 10~ 12 fj,' w . 
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where {. . represents coupling of rank one tensor oper- 
ators D and S^) to an operator of rank I. This coupling 
is a structure common to any PRCC term of i?l^^P. 
From the triangular condition, I — 0,1,2 are the al- 
lowed values, however, what values of I contribute de- 
pends on j v and j w . For example, I = 0, 1 contribute in 
the PNC 6 2 S 1/2 -> 7 2 S 1/2 transition of atomic Cs [J, 
where as only I — 2 contributes to the proposed PNC 
6 2 Si/2 — > 5 2 D 5 / 2 transition in Ba + [3]. 

The angular momentum diagram in Fig. [2jb), after 
evaluation, reduces to a 9j-symbol and free line part. 
Algebraically, the matrix element in the hyperfine states 
is 
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TABLE II. Component wise contribution from the coupled 
cluster terms for 6 2 S\/ 2 — > 7 2 Si/ 2 transition in Cs 
6 2 Si/2 - *• 5 2 D 3 / 2 transition in Ba 
transition in Ra + . 
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iM Cs 


3 


3 


-0.278 


4.198 


-0.005 


-1.533 
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4 


-0.317 


4.779 


-0.006 


-1.746 
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3 


0.764 


6.385 


-0.486 


-0.531 




3 


4 


0.657 


8.001 


-0.488 


-1.121 


136 Ba+ 


3 


2 


-2.676 


0.590 


-0.647 


-0.301 




2 
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2.566 


-0.723 


0.603 


-0.276 


139 Ba+ 
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-6.808 


1.729 


-1.620 


0.747 
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1.639 


-0.763 
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where D cS = ^{[xf^DS^T} 1 , is the effective dipolc 
operator in the hyperfine states. As seen from the an- 
gular momentum diagram, coupling of angular momenta 
in proper sequence is essential to obtain correct angular 
factors. However, the sequence is not manifest in the 
algebraic expression. 



(13) 



Results. — For the calculations reported in the letter, 
we use Gaussian type orbitals generated with V N_1 cen- 
tral potential. The iJlpj^, of Cs, Ba + and Ra + between 
various hyperfine states are given in Table. [TJ There is a 
close match between our MBPT results and results from 
similar works. 

There are changes when the transition amplitudes are 
calculated with PRCC. This can be attributed to the 
inclusion of higher order correlation effects. However, 
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it require a systematic series of calculations to examine 
the nature of the correlation effects from the higher or- 
der terms which are subsumed in the PRCC calculations. 
The results of 229 Ra + is a cause for concern, there is a 
large cancellation in the Fi = 3 — > Ff = 2 transition 
amplitude. However, for the other two transitions of the 
same ion, the transition amplitudes are higher than Ba + . 
In particular, the Ft = 2 — > Ff — 2 transition amplitude 
of 229 Ra + is the largest among all the values and this 
is in agreement with the previous results. For the neu- 
tral atom Cs, the PRCC results are larger than MBPT. 
This indicates, higher order correlation effects enahances 
SlpNc- R is opposite in Ba + and Ra+, the PRCC results 
are lower than MBPT and indicates higher correlation ef- 
fects have suppression effect. 

To examine the impact of electron correlation in bet- 
ter detail, consider the leading order (LO) and next to 
leading order (NLO) terms as listed in Table, [Il] In the 
PRCC calculations, as given in Eq. ( 11 ), for Cs these are 

D and DS^\ respectively. Here, the former repre- 
sents -ffpNc perturbed 7 2 Si/ 2 and has larger opposite 
parity mixing as it is energetically closer to odd parity 
states like 6 2 P\/2- The same is not true of 6 2 <Si/2, which 

is represented by DS^K 

In the case of Ba + the LO and NLO are 



Si D, respectively. Although, not shown in Table, 
similar pattern is observed in Ra + . The sequence is op- 
posite to Cs. Reason is, the transitions in these ions are 
of n 2 Si/2 — > n' 2 D 3 / 2 type and matrix elements of -ffp^c 
involving n' 2 D 3 /2 are negligible. Dominant contribution 
arises from the sp matrix elements, which are large. So, 
the term DS^, which represents ifpNc perturbation of 
n 2 Si/2 is the LO term of these ions. The contribution 

from D is, however, non-zero as n' 2 D 3 / 2 acquires 
opposite parity mixing through electron correlation ef- 
fects. It must be mentioned that, the Dirac-Fock contri- 
bution is the most dominant, however, in PRCC it is sub- 
sumed in the LO and NLO terms. The terms which are 
second order in cluster operators, in Eq. (11), are non- 



and 



II 



zero but small. For comparison, the two dominant con 

tributions from the second order term, DS^ 
it's hermitian conjugate, is given in the Table. |ll| 

2 



and 



We have also calculated the £Tp^c of the 6 2 Si 12 



5 2 L> 5/ 2 and 7 2 S 1 



/2 



6 -D5/2 transitions in Ba + and 



Ra + , respectively, and the results are given in Table. Ill 
The results from the PRCC are much larger than the 
MBPT results and this shows, without any ambiguity, 
electron correlation is the key to get meaningful results. 
This is on account of d 5 / 2 in the atomic n 2 D 5 / 2 states, 
which arc diffused and leads to larger electron correla- 
tions. 

Conclusions. — The PRCC theory we have developed 
incorporates electron correlation effects arising from a 
class of diagrams to all order with a nuclear spin- 
dependent interaction as a perturbation. It is a suitable 



TABLE III. Reduced matrix element, £1£nc> of tne 
6 2 Si/2 5 2 D 5 /2 and 7 2 Si/ 2 — > 6 2 D 5 / 2 transitions be- 
tween different hyperfine states Ba + and Ra + respectively. 

2 / 



The values listed are in units of ieag x 10 
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This work 




Other works 
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PRCC 




i3b Ba+ 
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0.003 


0.098 


0.227 


0.041 [19] 
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1 


0.002 


0.048 


0.127 




139 Ba+ 


2 


3 


0.003 


0.125 


0.235 


0.043 OH 


223 Ra+ 


3 


2 


-0.040 


-0.605 


1.262 


-0.526 HSU 


229 Ra+ 


2 


3 


-0.019 


-0.324 


0.616 


-0.256 QI] 



theory for precision calculations of atomic PNC arising 

rNSD 
Z PNC- 



from iJp^Q. With this method, it is possible to incorpo- 
rate electron correlation effects within the entire config- 
uration space obtained from a set of spin-orbitals. 
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